Chiral asymmetry of the Fermi surface in dense relativistic matter in a magnetic field 
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It is revealed that in the normal phase of dense relativistic matter in a magnetic field, there exists 
a contribution to the axial current associated with a relative shift of the longitudinal momenta in the 
dispersion relations of opposite chirality fermions. Unlike the topological contribution in the axial 
current at the lowest Landau level, recently discussed in the literature, the dynamical one appears 
only in interacting matter and affects the fermions in all Landau levels, including those around the 
Fermi surface. The induced axial current and the shift of the Fermi surfaces of the left-handed and 
right-handed fermions are expected to play an important role in transport and emission properties 
of matter in various types of compact stars as well as in heavy ion collisions. 
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Introduction. — At zero density and temperature, the 
structure of the ground state in relativistic chiral in- 
variant theories in a magnetic field is dictated by the 
magnetic catalysis phenomenon: the magnetic field is a 
strong catalyst of spontaneous chiral symmetry breaking 
[3, • The situation becomes much more complicated in 
the case of dense relativistic matter in a magnetic field 
Q. The main goal of this Rapid Communication is to 
reveal and describe some universal properties of such dy- 
namics. 

The recent studies [3] of similar dynamics in graphene 
in 2 + 1 dimensions have revealed several types of or- 
der parameters whose analogs have not been discussed 
in the context of relativistic models in 3 + 1 dimensions. 
(For earlier applications of the magnetic catalysis phe- 
nomenon in graphene, see Rcf. Q.) This motivated us 
to reexamine the properties of dense relativistic matter 
in an external magnetic field in 3+1 dimensions. As 
we show, an external magnetic field induces two quali- 
tatively different contributions to the net axial current 
in the normal phase of dense relativistic matter. The 
first one is a topological contribution due to the lowest 
Landau level (LLL) that was previously discussed in the 
literature [1, 0,(1, 0]. It exists even in free theories. (For 
related studies in hot quark-gluon plasma, see Ref. [lfl]-) 
In this Rapid Communication we find that there is also 
an additional dynamical contribution to the axial current 
that appears only in interacting matter. The crucial new 
point is that all Landau levels, including those around 
the Fermi surface, contribute to this interaction-driven 
contribution. Its origin is related to a relative shift of 
the longitudinal momenta in the dispersion relations of 



'Electronic address: gorbar@bitp.kiev.ua 

^Electronic address: vmiransk@uwo.ca On leave from Bogolyubov 
Institute for Theoretical Physics, 03680, Kiev, Ukraine, 
t Electronic address: igor,shovkovy@as u.edu| 



opposite chirality fermions. The amount of the shift is 
proportional to a coupling constant, the magnetic field, 
and the fermion chemical potential. Notably, it is almost 
independent of temperature. As we discuss below, this 
effect may have profound implications for the physics of 
compact stars as well as heavy ion collisions. 

Model. — In this paper, in order to illustrate the effect 
in the clearest way, a simple model will be utilized. We 
perform a study of a Nambu-Jona-Lasinio model with 
the Lagrangian density 

C = i> {iD v + fx S° v ) 7 > + ^ [{Uf + $i7V) a 

"(1) 

where 7 5 = Z7°7 1 7 2 7 3 and hq is the chemical potential. 
This model possesses the chiral U(1)l x U(1)r symmetry. 
The covariant derivative D v = d v — ieA v is taken in the 
Landau gauge, i.e., A v = xB8^ where B is the strength 
of the external magnetic field pointing in the z-direction. 

The structure of the full fermion propagator is given 
by a (3+l)-dimensional generalization of the ansatz used 
for the description of the quantum Hall effect dynamics 
in graphene [J], namely 

(«,«') = Uid t + m)7° - O ■ 7) - 7T 3 7 3 

+ i/r/V + A 7 3 7 5 - m 6 4 (u-u'), (2) 

where 717. = i(dk — ieA^) is the canonical momentum 
and u = (t,r). The above propagator contains several 
dynamical parameters that are absent in the tree level 
propagator, 
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(3) 

The physical meaning of the parameters m and \x is 
straightforward: m is the Dirac mass and \x is the full 
chemical potential in interacting theory. From the struc- 
ture of the inverse propagator in Eq. (J2J), it is clear that 
fi plays the role of the anomalous magnetic moment [in 



2 



graphenc, it also can be interpreted as a chemicalpoten- 
tial related to a conserved (pseudospin) current Q]. 

The meaning of the last parameter, A, is more subtle. 
In the context of graphene, A is a mass parameter that is 
odd both under time-reversal and parity transformations. 
In 2+1 dimensional models, this mass is responsible for 
inducing the Chern-Simons term in the effective action 
for gauge fields [ll[. In 3+1 dimensions, as suggested by 
Eq. |[5J), it is related to an induced axial current along 
the direction of the magnetic field, -07 3 7 5 -0. As will be 
shown below, A determines a shift of the longitudinal 
momenta in the dispersion relations of opposite chirality 
fermions in the chiral limit. We will call it a chiral shift 
parameter. Let us emphasize that, in the presence of 
an external magnetic field, the time-reversal and parity 
symmetries are broken and a state with the vanishing A 
is not protected by any symmetry. 

The parameters m, [i, A and jl are self-consistently 
determined from the gap equation, which takes the fol- 
lowing form in the mean field approximation: 

G- l {u,v!) = S- 1 (u 1 u')~iG int {G(u,u)~ 1 5 G(u,u) 1 5 
- tr[G(«,u)] +7 5 tr[ 7 5 G(u,w)]}(5 4 (w-M')- 

(4) 

While the parameters A and [i do not break the chiral 
U(1)l x U(1)r symmetry, nonzero values of m and jl 
break it down to U{\)l+r- In the mean field approxi- 
mation used here, we find that jl = in a self-consistent 
solution to the gap equation. While this fact simplifies 
the analysis, we note that jl may be non-vanishing in 
more refined approximations and in models with other 
types of interactions 0, [l2j ■ On the other hand, as sug- 
gested by the analysis in graphene, a nonzero jl should 
not change the main qualitative features of the phase 
with an induced A Q ■ 

The spectrum of fermionic quasiparticles is determined 
by thepoles in the full propagator ([2]) with jl = 0. As in 
Refs. 0, 3, expanding the propagator over the Landau 
levels, we arrive at the following dispersion relation: 



nonzero chemical potential. However, when fj,Q exceeds 
a certain critical value /x c , the chiral symmetry restora- 
tion and a new ground state are expected. As we show 
here, the new state is characterized by a no n- vanishing 
chiral shift parameter A and a nonzero axial current in 
the direction of the magnetic field. Since no symmetry 
of the theory is broken, this state describes the normal 
phase of matter that happens to have a rather rich chiral 
structure. 

The value of the Dirac mass mo in the vacuum state 
was calculated in Ref. In the weakly coupled regime, 
g = Gi nt A 2 /(47T 2 ) <C 1, the solution reads 
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where 1 = 1/ y/\eB\ is the magnetic length and A is the 
ultraviolet cutoff in the model at hand. This zero tem- 
perature solution exists for fi a < m . In this solution, 
the full chemical potential fi = /j,q. 

Our analysis shows that, in addition to the solution 
with a nonzero Dirac mass m, the gap equation also al- 
lows a solution with m = and a nonzero chiral shift 
parameter A, 



A = gfj, 



eB 



A 2 [l + 2ag]+g\eB\ 



(7) 



where a is a dimensionlcss constant of order 1, which is 
determined by the rcgularization scheme used in the anal- 
ysis (see below). Interestingly, the temperature depen- 
dence of A comes only through the chemical potential, 
which has a weak temperature dependence when T <C fi 
[l3T |. At T = 0, the chemical potential satisfies the fol- 
lowing equation: 



A* = Mo 
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V^ 2 - 2n\eB\6 (^ 2 - 2n\eB\) (8) 
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k 2 + a A sign(eB) 
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where n labels the Landau levels, a = ±1, and ks is the 
momentum along the direction of the magnetic field. In 
the chiral limit with m = 0, the states with the quan- 
tum number a = ±1 have a shifted longitudinal momen- 
tum, fc 3 — > fc 3 ± Asign(eB). As will be shown below, 
in this limit, the two different values of a correspond to 
fermions with opposite chiralities. As a result, the Fermi 
surfaces for left-handed and right-handed fermions be- 
come shifted. 

Results. — In accordance with the magnetic catalysis 
scenario [l|, Q , the ground state in the NJL model at van- 
ishing fiQ is characterized by a nonzero Dirac mass mo 
that breaks the chiral J7(1)l x U(1)r symmetry. Such 
a vacuum state can withstand a finite stress due to a 



The approximate solution to this equation is fi ~ (j,q up 
to power corrections in small g. 

The free energies of the two states are given by [l3| 
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respectively. In deriving the last expression, we used 
the approximate relations [i ~ [Iq and A ~ g^eB / K 2 . 
By comparing the free energies in Eqs. ([9]) and (fT0|) , we 
see that the ground state with a nonzero A becomes fa- 
vorable when fiQ > mo/y/2. This is analogous to the 
Clogston relation in superconductivity [Til ]. 
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The properties of the two types of quasiparticles corre- 
sponding to a = ±1 in the dispersion relation in Eq. ([5]) 
with m = are further clarified by the structure of the 
full propagator: 

oo 

G(u, u) = G- a V- + {GnT>- + G+V+) , (11) 

n=l 

where V± = \ [l ± i^ 1 ^ 2 sign(eB)] are the spin projec- 
tors. (For vl ^ u, the propagator will be presented else- 
where . ) The functions with n > are given by 

± _ i\eB\-y° f diodk 3 

Lu + /j±[k 3 - Asign(eB)] 
X I (w + n) 2 - 2n\eB\ - [fc 3 - A sign(eB)] 2 5 
uj + /i + [fc 3 + A sign(eff)] ,- 
(w + m) 2 - 2n|ftB| - [fc 3 + A sign(eS)] 2 5 J ' 

(12) 

where = | [l ± 7 5 sign(eB)] are the chirality projec- 
tors for a fixed sign of eB. As follows from this equa- 
tion, the quasiparticles of opposite chiralitics have dis- 
persion relations that differ from those in the free the- 
ory by the shift of their longitudinal momentum fc 3 — > 
fc 3 ± Asign(cB). This has profound implications for the 
physical properties of matter. 

The ground state with A ^ is characterized by a non- 
vanishing expectation value of the axial current density, 

(f 5 (u)) = -tr[ 7 VG M ]=^u-gU 

I rjl OO 

- ^AW v /2^p|,A) ! (13) 

n— 1 

where k(x, A) is a smooth cutoff function defined by the 
value of the cut-off energy A and a certain width of the 
region in which the value of the function drops from 1 
down to (i.e., k(x,A) ~ 1 for i « A and n(x,A) ~ 
for x ^ A). Taking this into account, we find that 

E^=i K{V 2n \ eB l A ) = aA 2 /\eB\ and a = 0(1). 

While the first two terms in the current density (jT3|) 
come from the LLL, the last term is due to the higher 
Landau levels, n > 1. Notably, all higher Landau levels 
below the (smeared) cut-off energy give (nearly) identical 
contributions to the induced axial current. The contri- 
bution proportional to /i is topological in nature and ap- 
pears even in the free theory p}. All other terms, which 
are proportional to A, are the result of interactions and 
have not been revealed in the literature before. 

We expect that the interaction-driven contributions 
can strongly modify transport and emission properties 
of dense relativistic matter in a magnetic field. Indeed, 
the corresponding induced axial currents are the conse- 
quence of a relative longitudinal flow of opposite chiral- 
ity quasiparticles, including those in higher Landau levels 
around the Fermi energy. This is in contrast to the role 



of the topological contribution that is exclusively due to 
the LLL, which is to a large degree quenched from the 
low-energy dynamics by the Pauli exclusion principle in 
many realistic cases. 

Discussion and Summary. — In this Rapid Communi- 
cation we found that, in accordance with the magnetic 
catalysis scenario, the vacuum state of relativistic matter 
in a magnetic field is characterized by a non-zero Dirac 
mass given by Eq. ^ at weak coupling. However, when 
the chemical potential exceeds a certain critical value, 
/i c ~ nio/V2, such a state is replaced by the normal 
ground state that is characterized by the following two 
properties: (i) the presence of an induced axial current 
along the magnetic field and (ii) the presence of the dy- 
namically generated chiral shift parameter A, which is a 
3 + 1 dimensional analog of the parity odd mass term in 
2 + 1 dimensions leading to the Chern-Simons term. We 
find that, in addition to the previously known topological 
term in the induced current, there are also interaction- 
driven contributions from the lowest as well as from the 
higher Landau levels. In fact, the newly found contri- 
butions arc directly related to a dynamically generated 
value of the chiral shift parameter A ~ g/xoeB/A 2 . This 
parameter quantifies the relative shift of the longitudinal 
momenta in the dispersion relations of opposite chirality 
quasiparticles. 

It might be appropriate to mention that, for instruc- 
tional purposes, our study here is simplified: we used an 
NJL-type local interaction and utilized the mean-field ap- 
proximation. These limitations may lead to the results 
that are not always quantitatively reliable, e.g., in the 
context of stellar matter. Nevertheless, it is expected 
that our results should remain qualitatively the same 
even when more realistic models are used. Indeed, the 
fact that the expression for A in Eq. ([7]) is linear in g 
in the lowest order indicates that the corresponding dy- 
namics is essentially perturbative. Apparently this is a 
general feature that should not depend on whether the 
interactions are short range, as in the NJL model, or long 
range, as in QCD or QED. In either case, a vanishing A 
is not protected by any symmetry. 

Another limiting assumption of this study is the exact 
chiral symmetry. However, most of the results are not 
modified much when the symmetry is at least approxi- 
mate, i.e., when the fermions have non-zero bare Dirac 
masses, but such masses are small compared to the value 
of the chemical potential For the applications in 

protoneutron stars suggested below, this approximation 
is justified, but the general study will be of interest, e.g., 
in relation to the electron plasma in white dwarfs. 

In the future, it will be also of interest to address a 
possible interference of the dynamics responsible for the 
generation of the chiral shift parameter with color super- 
conductivity in quark matter [l5| . Here we just mention 
that the normal ground state with a nonzero A seems to 
be the only possibility at higher temperatures, which are 
of main interest for us here. 

We expect that the generation of the chiral shift pa- 
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ramcter may affect physicaf properties of the quark mat- 
ter in quark and/or hybrid stars, the electron gas in neu- 
tron stars, and possibly even the electron gas in white 
dwarfs. The corresponding fcrmionic systems are degen- 
erate (T -C fJ>) and the results of this Rapid Communica- 
tion apply directly. It is appropriate to mention, however, 
that the chiral shift parameter does not vanish even in 
the non-degenerate limit (T> /i), although the analysis 
of the dynamics becomes more involved Therefore, 
the generation of a non-zero A can also affect the chiral 
magnetic effect in heavy ion collisions [l(| . 

One of the consequences of the phenomenon discussed 
in this Rapid Communication is the possibilit y o f a qual- 
itatively new mechanism for the pulsar kicks [l6[ . In the 
presence of a magnetic field, almost any type of relativis- 
tic matter inside a protoncutron star should develop ax- 
ial currents as in Eq. (|13[) . The main carriers of such 
currents are the electrons in the nuclear matter, and 
the quarks together with the electrons in the deconfincd 
quark matter. Since the induced currents and the chiral 
shift parameter have only a weak temperature depen- 
dence (assuming T <^ fi), this phenomenon may provide 
a robust anisotropic medium even at the earliest stages 
of the protoncutron star. This is important because 
moderately hot matter with 10 McV < T < 50 MeV, 
present during the first few tens of seconds of the pro- 
toneutron star evolution [l7| . may have a large enough 
amount of the thermal energy to power the strongest 
(with v > 1000 km/s) pulsar kicks observed [l||. In con- 
trast, the constraints of the energy conservation make 



it hard, if not impossible, to explain such kicks if the 
interior matter is cold (T < 1 MeV). The common diffi- 
culty of using a hot matter, however, is the very efficient 
thermal isotropization that washes out a non-isotropic 
distribution of neutrinos produced by almost any mech- 
anism [l8l . [l|| . In the mechanism proposed in this Rapid 
Communication, however, the asymmetric distribution of 
the neutrinos develops as a result of their multiple clastic 
scattering on the left-handed electrons or quarks, flowing 
in the whole bulk of the stellar matter in one direction 
along the magnetic field. 

In passing, let us mention that the robustness of the 
axial currents in hot magnetized matter may also provide 
an additional neutrino push to facilitate successful super- 
nova explosions as suggested in Ref. [2(| ■ Further details 
of this mechanism will be discussed elsewhere [l3[ . 
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